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We explore equilibrium solutions of non-topological solitons in a general class of scalar field 
theories which include global U(f) symmetry. We find new types of solutions, tube-shaped and 
crust-shaped objects, and investigate their stability. Like Q-balls, the new solitons can exist in 
supersymmetric extensions of the Standard Model, which may responsible for baryon asymmetry and 
dark matter. Therefore, observational signals of the new solitons would give us more informations 
on the early universe and supersymmetric theories. 



PACS numbers: 03.75.Lm, 11.27.+d 



I. INTRODUCTION 

In apioneering work by Friedberg, Lee and Sirlin in 
1976 1 , non-topological solitons were introduced in a 
model with a U(l)-symmetric complex scalar field cou- 
pled to a real scalar field. In contrast with topological 
defects, they are stabilized by a global U(l) charge, and 
their energy density is localized in a finite space region 
without gauge fields. In 1985 Coleman showed such soli- 
tons exist in a simpler model with an S0(2) (viz. Ufl)) 
symmetric scalar field only, and called them Q-balls 0- 

Q-balls have attracted much attention in particle cos- 
mology since Kusenko pointed out that they can exist 
in all supersymmetric (SUSY) extensions of the Stan- 
dard Model '3]. Specifically, Q-balls can be produced 
efficiently in the Affleck-Dine mechanism [4i] and could 
be responsible for baryon asymmetry Q and dark mat- 
ter . Q-balls can also influence the fate of neutron stars 
Based on these motivations, stability of Q-balls has 
been intensively studied [8l-[Tl|. 

Observational signatures of SUSY Q-balls has been 
studied [l^ . and their mass and flux were constrained 
by experimental data of the searches for magnetic 
monopoles and heavy cosmic rays [isj . Currently direct 
searches for neutral Q-balls and for electrically charged 
Q-balls are in progress in Super-Kamiokande II [l3| and 
in the SLIM Experiment 15| , respectively. Furthermore, 
it has been shown that gravitational waves are emitted 
during Q-ball formation and could be detected by next- 
generation gravitational detectors pLS], 

In spite of increasing concern about non-topological 
solitons in SUSY, other equilibrium solutions has not 
been studied so much, while topological defects have sev- 
eral types according to the symmetry. In this paper we 
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address a fundamental question: are their any other non- 
topological solitons in different forms? If new solitons 
exist in the theories which allow for Q-balls, their obser- 
vational signals would give us new informations on SUSY 
Q-ball models. Here we re-analyze scalar field theories 
which include global U(l) symmetry. 



II. EQUILIBRIUM SOLUTIONS 

We begin with a review of Q-ball solutions. Consider 
an S0(2)-symmetric scalar field 4> = {.4'ii4>2), whose ac- 
tion is given by 



S 



\ a=l 



Due to the symmetry there is a conserved charge. 



= d^xlcj) 



dt 



dt 



(1) 



(2) 



Assuming spherical symmetry and homogeneous phase 
rotation, 



4> — (/) (r) {cos ujt, sin ujt), 
one has a field equation, 
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dV 



dr"^ ~^ r dr ^ ^ dd) 
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This is equivalent to the field equation for a single static 
scalar field with a potential = V — uj^ (j)^ / 2 . 

Equilibrium solutions 4'{r) with a boundary condition 



dr 



(r = 0) = 0, </.(r ^ cx)) = 0, 



(5) 



exist if min(K;) < K;(0) and (iVL;/d<^(0) > 0. This con- 
dition is rewritten as 
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FIG. 1: Interpretation of (a) Q-ball solutions and (b) new soli- 
ton solutions by analogy with a particle motion in Newtonian 
mechanics. 



If one regards the radius r as 'time' and the scalar 
amplitude (f)(r) as 'the position of a particle', one can 
understand Q-ball solutions in words of Newtonian me- 
chanics, as shown in Fig. [TJa). Equation (jH) describes 
a one-dimensional motion of a particle under the con- 
served force due to the potential — and the 'time'- 
dependent friction —{2/r)d(j)/dr. If one chooses the 'ini- 
tial position' 0(0) appropriately, the static particle be- 
gins to roll down the potential slope, climbs up and ap- 
proaches the origin over infinite time. 

To demonstrate numerical solutions later, we adopt a 
simple model, 

V = ^m^ct)'^ - ficl)^ + A0^ with m^, ^, A > 0. (7) 



and rescale the quantities as 



m, uj — uj. 



A 

VA ^J■ 

Then, the existing condition (jG)) becomes 
0<.^<i, e- 



(8) 



(9) 



string-like configuration, 

= (f){R) {cos {nip + ujt),sin{mp + ujt)), (10) 

where n is nonnegative integer and {R,ip, z) is the cylin- 
drical coordinate system. The field equation becomes 



+ 



1 



dV 



(11) 



If n = 0, the field equation is the same as ^ except for 
a numerical coefficient. Therefore, Q-ball like solutions 
of (j){R) exist. If n > 1, there is no regular solution 
which satisfies 0(0) 7^ 0. However, if we adopt a different 
boundary condition. 



(p{R = 0) = 4>{R^oo) =0, 



(12) 



there is a new type of regular solutions. We introduce an 
auxiliary variable which is defined by 0(i?) ~ R^ip{R), 
Then, Ea. ([TT|) becomes 



d'^ip 2n + ldtp 
dR^ ^ R dR 



(13) 



If we choose ■0(0) appropriately, we obtain a solution 
?/'(i?) which is expressed in the Maclaurin series without 
odd powers in the neighborhood of i? = 0. In terms of 
the original variable 0(i?) , the nth differential coefficient 
0(")(O) — ^0(0) should be determined by the shooting 
method, while any lower derivative vanishes at i? = 0. 
We plot some solutions in Fig. [2lja). 

We can illustrate existence of the new solutions with 
^ 1 by analogy with a particle motion in Newtonian 
mechanics, as shown in Fig. [ijb). Equation (fTTj) de- 
scribes a one-dimensional motion of the particle under 
the conserved force due to the potential — K;(0) and two 
non-conserved forces, the friction —{l/R)d(j)/dR and the 
repulsive force n^0^/i?^. If n = 1, by choosing the 'ini- 
tial velocity' d(j)/dR{0) appropriately, the particle goes 
down and up the slope, and at some point = 0max it 
turns back and approaches the origin over infinite time. 
If n > 2, d(p/dR{0) vanishes; instead, the nth derivative 
0^"^(O) gently pushes the particle at = 0. Therefore, 
with the appropriate choice of 0'"' (0), the particle moves 
along a similar trajectory to that of n = 1. This argu- 
ment also indicates that the existence condition of the 
new soliton solutions are the same as that of Q-balls, © 
or Solutions with the same behavior as the n = 1 
solutions were obtained by Kim et aL[l3], who studied 
the S0(3)-symmetric scalar field without Q-charge. 

Nonzero components of the energy-momentum tensor 
are given by 



dR 



(1) Q-tubes. For the same S0(2) model, we suppose a 
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FIG. 2: Examples of Q-tube solutions, (a) shows </>(?-) for four 
cases, (b) and (c) show each term of the energy-momentum 
tensor for two of the solutions, where we put = 0.6. 



tion as follows. If we define the gravitational potential 
$ as $ = gtt + 1, and take it into account up to its first 
order, the Einstein equations yeild the extended Poisson 
equation. 



where T? is the trace of the spatial components. Except 
for n = cases, the gravitational source in (fT5|) vanishes 
in the center in the x-y plane. Furthermore, this field 
configuration has planar symmetry in the z direction. 
Therefore, we collectively call the new solitons Q-tubes. 

Because S0(2) is a subgroup of SO(iV > 3) or SU(iV > 
2), Q-tubes with n = as well as Q-balls can appear in 
any SO(A^) or SU(iV) theory. By contrast, we suspect 
that Q-tubes with n > 1 can appear only in U(l) the- 
ory, since their string-like configuration is topologically 
unstable in theories with SO{N > 3) or SU{N > 2). 
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The nonzero component indicates that the solutions 
possess angular-momentum. Each term of the energy- 
momentum tensor is presented in Fig. [2Ub)(c). Although 
the present analysis does not include gravity, we can es- 
timate gravitational effects in the weak field approxima- 



FIG. 3: Examples of Q-crust solutions, (a) shows <^(r) for two 
cases, (b) shows each term of the energy-momentum tensor 
for the solution with = 0.49. 



(2) Q-crusts. Next, we consider an S0(3)xU(l)- 
symmetric scalar field 4> = 02, ^3), whose action 

is given by 



5 = I d'^x 



A ^<f>a(l)a 
\ a=l 

(16) 
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Assuming 

(p — e*'^*(^(r)(cos(/3sin6',sin93sin6',cos0), (17) 
we obtain the field equation, 



(P<i) 2 d<i) 2(f) 



dV 



dr"^ ~^ r dr j.2 ^ J-i ( 



If we adopt a boundary condition, 

0(r = 0) = if){r oo) ~ 0, 



(19) 



we find a regular solution (f>{r), like a Q-tube solution 
with n = 1. We plot some solutions in Fig. |3Ia). 

Nonzero components of the energy-momentum tensor 
are given by 



2 V rfr 



1 (dcpy 02 



(20) 



In this case T^^ vanishes, which means that the solutions 
possess no angular- momentum. Each term of the energy- 
momentum tensor is also shown by Fig. [3l^b). The ex- 
tended Poisson equation is given by 

= AttG{-TI + T^) = %TTG{uj^(t>'^ - V). (21) 

We find that the kinetic term uP'ip'^ , which is responsi- 
ble for Q-charge, is dominant but vanishes in the center. 
Therefore, we call the solutions Q-crusts. 



III. PROPERTIES OF THE SOLUTIONS 

Stability of Q-balls has been discussed essentially by 
energetics. The total energy of the system is defined by 



1 



1 



d-'xl-Lo'cf)' + -{d^(t)Y +V 



(22) 



where di denotes a spatial derivative. For a fixed model 
V{(j)) and the phase assumption ([3]), there remains a free 
parameter, w or Q; accordingly, there is a family of equi- 
librium solutions. For such a family, if E increases as a 
function of Q but dE/dQ decreases, energetics prohibits 
one Q-ball splitting into two under fixed Q. In this case 
we can understand that Q-balls are stable under the as- 
sumption ([3]). In connection with this argument, Paccetti 
Correia and Schmidt showdeed a useful theorem that sta- 
bility is determined by the sign of {uj/Q)dQ/duj 0. 

Sakai and Sasaki [lOj proposed a simple method of an- 
alyzing stability using catastrophe theory [l^ as follows. 
Catastrophe theory reveals stability of a mechanical sys- 
tem completely once behavior variable(s), control param- 
eter (s) and a potential are given. Therefore, an essential 



point is to choose those variables in the Q-ball system ap- 
propriately. For a given potential V{(f>) and charge Q, we 
consider virtual displacement 5(j){r) near the equilibrium 
solution 4>{r). If we redefine lo by 



oj = q/ f <j>l{x)d^x, 



(23) 



the domain of definition of uj is extended to off- 
equilibrium configurations. Using this w, we can repre- 
sent a continuous deformation by a one-parameter fam- 
ily of displacement functions, 5(j)ui{r). Then the energy 
([22|) is regarded as a function of lo: E{uj) = E[(j)uj]. Be- 
cause dE/doj = {5E / 54>i^)d(j)i^ / dio = when (f)^ is an 
equilibrium solution, u may be regarded as a behavior 
variable and E as the potential. The charge Q is given 
by hand, or physically, it is determined by initial condi- 
tions; therefore, it should be regarded as a control param- 
eter. Catastrophe theory tells us that stability changes 
at dQ/dE = or dQ/duj — 0, which are consistent with 
the above arguments. 

Properties of Q-balls in the model ([7]) were elucidated 
as follows [lo| . 

• > 1/2: y(0) is the absolute minimum. There is 
no bounds on Q, and all equilibrium solutions are 
stable. 

• rh? < 1/2: 1^(0) is a local minimum but the abso- 
lute minimum is located at 7^ 0. For each m^, 
there is a maximum charge, Quia.x, above which 
equilibrium solutions do not exist. For Q < Qmax, 
stable and unstable solutions coexists. 

It turns out that properties of Q-tubes and Q-crusts also 
depend on whether > 1/2 or < 1/2. Therefore, 
in the following, we show numerical results for rh^ = 0.6 
and = 0.3 as typical examples. 

(1) Q-tubes. In this case, E and Q diverge because 
they are infinitely long. Therefore, we define the energy 
and charge per unit length, respectively, as 



2n 



q ~ 27rw 



^dR. 



2i?2 



V 



(24) 



In accordance with the normalization (|8]), we rescale the 
energy / charge variables as 



A2 xi -. \i 

e=— , q= — g, E= — E, Q^XQ. 
M A* M 



(25) 



We show the q-e relations for = 0.6 and for = 
0.3 in Fig. From a viewpoint of these relations, basic 
properties for Q-tubes are the same as those for Q-balls 
as described above. They are summarized as follows. 

• fh? > 1/2: There is no bounds on q, and all equi- 
librium solutions are stable. For fixed q, Q-tubes 
with lower n are energetically more stable. 
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FIG. 4: Existence domain of Q-tube solutions in q-S space FIG. 5: Existence domain of Q-crust solutions in Q-E space 

and their stability. We put = 0.6 and = 0.3 in (a) and their stability. We put = 0.6 and = 0.3 in (a) and 

and in (b), respectively. The solid and dashed lines represent in (b), respectively, 
stable and unstable solutions, respectively. 



• < 1/2: For each m^, there is a maximum 
charge, gmax- For q < qmax, stable and unstable so- 
lutions coexists for fixed q. Interestingly, Q-tubes 
with higher n can have larger q. 

(2) Q-crusts. Similarly, we show the Q-E relations for 
w? — 0.6 and for ttt? — 0.3 in Fig. [5l For reference, we 
show the results for Q-balls, too. (Note that Q-balls can 
appear not only in U(l) theories but also in a wide class 
of theories which include global U(I) symmetry.) What 
we find for Q-crusts is summarized as follows. 

• > 1/2: There is no bounds on Q, and all equi- 
librium solutions are stable. For fixed Q, Q-balls 



are energetically more stable than Q-crusts, as ex- 
pected. 

• < 1/2: For each m^, there is a maximum 
charge, Qmax- For Q < Qmax, stable and unsta- 
ble solutions coexists. Interestingly, Q-crusts can 
have larger Q than Q-balls. 



We should note that the above results are based on the 
phase assumption ([TOl) or ([17]). Stability against pertur- 
bations on the phase configuration ((T(I| or p7)) cannot 
be revealed by the present energetic analysis and should 
be studied by dynamical analysis. 
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IV. SUMMARY AND DISCUSSIONS 

We explore equilibrium solutions of non-topological 
solitons in a general class of scalar field theories which 
include global U(l) symmetry, and then we find new 
types of solitons: Q-tubes in U(l) theories and Q-crusts 
in S0(3)xU(l) theories. Because only n = Q-tubes 
have homogeneous phase like Q-balls, they can appear in 
a wide class of theories which include U (1) symmetry. 
Except for n = Q-tubes, there is a dip in kinetic energy 
in the center. In contrast with cosmic global strings or 
global monopoles, their gravitational mass can be finite 
without gauge fields. 

We also investigate stability of equilibrium solutions 
for the model ([7]) under the phase assumption ((TU]) or 
P?)l by calculating their charge and energy (or those per 
unit length). The charge-energy relations indicate that, 
if V^(0) is the absolute minimum, there is no bound on 
charge, and all solutions are stable. On the other hand, 
if V{0) is a local minimum but the absolute minimum 
is located at 7^ 0, there is a maximum charge, above 
which equilibrium solutions do not exist. For fixed charge 
below the maximum, stable and unstable solutions coex- 
ists . It is interesting that Q-tubes with higher winding 
number can have larger charge density and that Q-crusts 
can have larger charger charge than Q-balls. Stability 
against these phase configurations is beyond the present 
energetic analysis and should be studied by dynamical 
analysis, which is our next subject. 



Unlike Q-crusts, our Q-tubes solutions are infinitely- 
long and unrealistic in themselves. Nevertheless, Q-tubes 
are the more interesting because they can appear in the 
minimal (i.e., U(l)) models and several researchers have 
already performed numerical simulations in those models 
[T9| . Those simulations showed that filament structure 
appears just before Q-ball formation and maintain its 
shape for a certain time. We conjecture that such a fila- 
ment structure is semi Q-tubes. Furthermore, according 
to recent simulations of the collision of two Q-balls [l^l , 
apparent two rings are formed. We suspect that they 
are loop Q-tubes. Further fundamental investigations of 
Q-tubes and the other solitons together with advanced 
simulations of Q-ball formation may confirm the above 
conjectures and elucidate observational consequences of 
Q-ball formation in SUSY. It may be also interesting to 
explore analogous new solitons in non-relativistic atomic 
Bose-Einstein condensates pTj . 
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